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Abstract: State of the art X-ray imaging sensors comprise a trade-off between the achievable 
efficiency and the spatial resolution. To overcome such limitations, the use of structured and scin- 
tillator filled aluminum oxide (AlOx) matrices has been investigated. We used Monte-Carlo (MC) 
X-ray simulations to determine the X-ray imaging quality of these AlOx matrices. Important fac- 
tors which influence the behavior of the matrices are: filling factor (surface ratio between channels 
and 'closed 4 AlOx), channel diameter, aspect ratio, filling material etc. Therefore we modeled the 
porous AlOx matrix in several different ways with the MC X-ray simulation tool ROSI and 
evaluated its properties to investigate the achievable performance at different X-ray spectra, with 
different filling materials (i.e. scintillators) and varying channel height and pixel readout. In this 
paper we focus on the quantum efficiency, the spatial resolution and image homogeneity. 
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1. Goals and motivation 

State of the art X-ray imaging sensors comprise a trade-off between the achievable efficiency and 
the spatial resolution. For instance a thick scintillator yields a high efficiency but suffers from a de- 
crease in spatial resolution due to (optical) light spreading in the scintillation volume. For thinner 
scintillators the opposite is observed. In general an increase in spatial resolution leads to a decrease 
in efficiency and vice versa. The goal of the Fraunhofer-internal project WISA Honeris is to over- 
come such limitations and to provide 'ready for the market' X-ray imaging sensors which offer 
both high efficiency and high spatial resolution. This is achieved by filling structures (hexagonally 
arranged cylindrical channels) of porous self-organized AlOx matrices with appropriate scintillator 
materials. Due to its channel-like structures which act as light guides, the high spatial resolution 
is maintained with increasing scintillator thickness. This permits the fabrication of very thick and 
therefore highly efficient scintillator matrices without losing spatial resolution. We will not go into 
further detail with the fabrication of the AlOx matrices, since this is not the focus of this paper. 
More information can be found in [0]. Instead we used Monte-Carlo simulations to investigate the 
effects of varying matrix parameters and different scintillation materials regarding spatial resolu- 
tion and quantum efficiency. The goal of this simulation study is to generate a priori knowledge 
about the behavior of the scintillator matrix and to examine various questions which are difficult to 
answer in real experiments. This should help in the fabrication and optimization of the matrices. 

2. Monte-Carlo model 

We constructed a model of the matrix structure in ROSI with adjustable parameters like channel 
diameter d c , channel height h c , channel pitch (center to center) p c and scintillation material. For 
the following simulations a d c of 1.6 jam and a p c of 2 jam were used. To reduce the simulation time 
we approximated the matrix structure by a simplified model consisting of a massive scintillator 
medium divided by diagonal and horizontal strips of certain thickness (see figure |T|). The filling 
factor was about 36 %. The irradiation was simulated with a tungsten X-ray tube at varying tube 
voltages U t . 
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Figure 1: Simplified simulation model. 

3. Quantum Efficiency 

The simulations regarding quantum efficiency (QE) were carried out with different scintillation 
materials, namely Gd 2 2 S (GOS), CdW0 4 (CWO), PbW0 4 (PWO) and Lu 2 Si0 5 (LSO). To de- 
termine the absorption properties of the matrix structure consisting of the scintillator channels and 
the surrounding AlOx volume, an area of 4 x 4 pixels with a size of 4 um 2 each and h c of 100 um 
and 500 um was irradiated with both U t = 35 kV and U t = 160 kV. Based on the deposited energy 
in the scintillator channels, we calculated the respective efficiency (excluding loss due to optical 
effects and limited readout filling factor, which were not simulated). The absorption probability 
curve for the homogeneous materials with h c = 100 um is shown in figure 0. Although the heavy 
materials PWO and CWO show a better mean probability over the whole energy range, GOS still 
reaches a good overall efficiency due to its (Gadolinium) K-edge at about 50 keV, which is close 
to the K-a emission line of the applied tungsten X-ray spectra. The absorption properties of the 
whole scintillator matrix are shown in table |T|. As expected, the U t = 35 kV case showed nearly 
similar results for all tested scintillator materials. With a QE near 100 % for all scintillators at the 
center of mass energy of the 35 kV spectrum at about 15 keV (see figure @), the efficiency is mostly 
determined by the AlOx volume which is identical in all cases. The deviation between the mate- 
rials is below 4 %. Simulations with U t = 35 kV at h c = 500 um were skipped because of similar 
results to be expected. With U t = 160 kV the differences were more significant. PWO performed 
best with 7.8 % (23.9 %) efficiency while GOS reached 6.9 % (20.2 %) for h c = 100 um (500 um). 
The results for CWO were quite similar with 6.9 % (22.6 %). LSO showed the worst QE with only 
5.2% (18%). 
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Table 1: Quantum efficiency for the AlOx matrix filled with different scintillator materials, as a 
function of h c and U t . The simulations for h c = 500 um at U t = 35 kV were skipped because of the 
negligible effect of the scintillator thickness on the efficiency at these low energies. The numbers 
are given in percent. 
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4. Image homogeneity 

To examine the influence of the number of channels per readout pixel on the image homogeneity, 
the hexagonal structure was read out by an orthogonal pixel matrix with varying pixel sizes ranging 
from 1 |um to 30 jam. The scintillator matrix was irradiated with U t = 40 kV. The superposition 
of two periodical grids, in this simulation setup, the scintillator matrix and the pixelated readout 
(sensor) leads to an intensity modulation as seen in figure [3]. The modulation is about 15 %. This 
is due to aliasing errors (Moire-pattern). The frequency of the modulation varies with the number 
of channels per readout pixel. The more channels per pixel the lower the modulation frequency. 
One phenomenon, which was observed during these simulations, was an intensity gradient from the 
center towards the outside of the scintillator matrix. This effect becomes more significant with an 
increasing number of channels per readout pixel. To take a closer look at this effect, the scintillator 
matrix was simulated with a fixed readout pixel matrix of 20 jam pixel pitch and h c of 100 jam and 
500 jam, irradiated with U t of 40, 80, 120 and 160 kV. We used a virtual Focus-Detector-Distance 
(FDD) of 100 cm and irradiated one half of the scintillator matrix of 4.8 cm edge length with a 
cone beam. This results in X-rays impinging the matrix at an angle 6 between 0° and 1.4° with 
respect to the normal of the matrix surface (see figure f|). 




Figure 3: Intensity modulation and Figure 4: Irradiation setup for the simu- 

inhomogeneities with different pixel lations regarding spatial resolution. Not 

sizes. true to scale. 
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Unlike a homogeneous scintillator medium where the intensity decreases towards the edges 
of the medium, a significant drop of intensity can be seen in the center of the simulated matrix. 
This intensity gradient depends on U t and h c as can be seen in figure |3] and The gradient results 
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Figure 5: Intensity gradient for 
h c = 1 00 jam. 
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Figure 6: Intensity gradient for 
h c = 500 jam. 



from X-ray photons impinging the matrix at an angle. More precisely, it depends on the integral 
scintillator thickness, which an X-ray photon penetrates on its trajectory through the matrix until 
absorption or leaving of the structure. In the simulated setup (fixed channel diameter and fixed pixel 
pitch) this in turn depends on h c and the X-ray photon energy. With h c = 100 jam and a low mean 
photon energy like with U t = 40 kV, the mean penetration length is rather small, which means the 
absorption happens mainly at the surface of the scintillator channel. For simplification one could 
think of the surrounding AlOx material as zero absorbing and the scintillation channels as fully 
absorbing. Then photons at a very small angle have a higher chance to pass the scintillator matrix 
without hitting a scintillator channel (see figure This results in a significant drop of absorbed 
photons, i.e. less intensity in the area near the center. Increasing h c reduces this angle where no 
absorption happens. This results in a steeper curve as can be seen in figure With higher photon 
energies the assumption that absorption mainly takes place at the surfaces is no longer true. The 
probability to be transmitted through scintillator channels increases. Furthermore photons are more 
likely to interact in adjacent channels if the incident angle is high enough. The integral scintillator 
thickness gets more important. For h c = 100 jam this effect results in a less significant drop of 
intensity near the center with increasing energy. The h c = 500 jam case shows an inverse effect. 
This means there must be another effect which affects the gradient. Although not investigated in 
detail we assume, that this has to do with beam hardening. The AlOx parts with h c = 500 jam can 
no longer be seen as zero absorbing. 

Imagine two photons p\ and p2. p\ impinges the scintillator matrix at such a small angle 
that it will go through the matrix without hitting a scintillator channel on its trajectory for both 
h c = 100 jam and 500 jam. On the contrary p2 impinges the matrix at a slightly larger angle, so that 
it only hits a scintillator channel if h c is at least 500 jam. In the 100 \im case the energy E of the 
particles p\, p2 is not important. They either get absorbed in the AlOx volume or are transmitted 
through the matrix. In both cases no signal is generated. The 500 jam case however shows an 
energy dependence for p2. For a small E p 2, P2 is more likely to get absorbed in the AlOx volume 
because of the long distance to reach the scintillator channel. In contrary if E p 2 is high, p2 will 
probably reach the channel and deposit some energy. That means the difference between p\ and 
P2 increases with increasing photon energy as can be seen in figure An X-ray spectrum with a 
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higher mean photon energy results in a steeper curve in the small angle area around the center. 
5. Basic spatial resolution 

We used the basic spatial resolution (BSR) method to study the spatial resolution as a function 
of the scintillation material, U t and h c . It is part of a well established detector characterization 
standard [BD- Compared to an MTF calculation, the BSR method is a more practical approach 
for measuring the visual spatial resolution. It utilizes a plate of acrylic glass containing pairwise 
arranged platinum and tungsten wires of decreasing diameter. It is placed directly in front of a de- 
tector. From the gray value profile perpendicular to the wire length, one can calculate the contrast 
of the wires in relation to the area between the wires. The BSR value is defined as the (interpo- 
lated) wire diameter with exactly 20 % contrast. Considering the very high spatial resolution to be 
expected, we modeled a modified double wire phantom with thinner wire diameters than specified 
in the standard, ranging from 8 jam to 1 jam. The simulated scintillation materials were GOS with 
the highest light yield of all considered materials 0], [0] and PWO for the highest QE (see table |T|). 
Two different irradiation setups were taken into account. In addition to the cone beam irradiation 
with 9 between —1.4° and +1.4° we also investigated a parallel beam geometry, where every wire 
pair is irradiated under the same angle = 0°. The cone beam setup corresponds to a realistic setup 
with an FDD of 100 cm and a scintillator matrix edge length of 4.8 cm. Each case was simulated 
with U t = 35 kV, 160 kV and with h c = 100 jam, 500 jam. The simulated image contrast for each 
wire pair can be seen in figure 0, the corresponding calculated BSR values in table @. With an X- 
ray spectrum of U t = 35 kV, corresponding to a small penetration length, we observed only small 
differences between the two scintillators. At parallel beam irradiation both materials show nearly 
the same spatial resolution with both h c = 100 jam and 500 jam. However, for cone beam irradia- 
tion the spatial resolution is significantly degraded (up to 17 % compared to the parallel case) by 
increasing h c from 100 jam to 500 jam because of the higher probability of X-rays impinging the 
matrix structure at an angle, to hit a scintillator channel. The difference between both materials is 
below 6 % due to the low Compton scattering cross section in this X-ray energy range. 

An X-ray spectrum of U t = 160 kV, corresponding to a higher penetration length, leads to 
an overall lower spatial resolution than the 35 kV case because of the generally higher Compton 
scattering cross section. Furthermore the higher X-ray energy increases the differences between 
the two tested materials (up to 9 %). At parallel beam irradiation only a small degradation can be 
seen when increasing h c from 100 jam to 500 jam. At cone beam irradiation the spatial resolution 
depends on both material and channel height. 

The material dependency is caused by the different X-ray interaction cross sections. As can be 
seen from figure 0, at U t = 35 kV PWO reaches slightly better contrasts than GOS due to its higher 
absorption capability while at U t = 160 kV GOS clearly surpasses PWO. We assume that this is 
directly related to the high energy X-ray fluorescences of Pb (75 keV), which occur at energies 
above 88.3 keV and are isotropically emitted. Compared to the Gd fluorescences of 43 keV, they 
have a much higher penetration length and therefore a higher impact on the degradation of the 
spatial resolution. 

The h c dependency is caused by the higher probability to absorb a scattered X-ray photon for higher 
channels. 
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Figure 7: Image contrast based on BSR method for GOS and PWO with (a) parallel beam and 
U t = 35 kV, (b) parallel beam and U t = 160 kV, (c) cone beam and U t = 35 kV, (d) cone beam and 
U t = 160 kV. The local minima at the 4 jam wire pair in figure a and b result from aliasing errors 
caused by the pixel readout. 
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Table 2: BSR values for the scintillator materials GOS and PWO, as a function of irradiation 
setup, h c and U t . The numbers are given in um. 

6. Conclusion and outlook 

The goal of this simulation study was to get insight in the behavior of structured scintillators with- 
out the need to fabricate multiple test samples. Therefore we used a simplified simulation model 
and studied the efficiency and spatial resolution with different materials and geometries. The sim- 
ulations clearly showed the achievable performance of this technique but also gave insight in the 
forthcoming problems and potential solutions. GOS turned out to be the most suitable of the tested 
materials due to its high spatial resolution, light yield and decent QE. Heavy materials like PWO 
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could not compensate their far inferior light yield with their high QE which is mainly determined 
through the dominating AlOx parts. In addition they did not show a better spatial resolution though 
their lower scattering contribution as shown in the BSR simulations. 

Further research has to be done on the topic of aliasing errors, which are to be expected because 
of the special channel-like geometry including methods to minimize these. In the meantime first 
prototypes based on the simulation results have been produced and are experimentally investigated 
and compared to the simulations. This will be subject to an upcoming publication. 
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